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Modeling the two-point correlation 
of the vector stream function 

By M. Oberlack 1 , M. M. Rogers 2 AND W . C. Reynolds 3 


A new model for the two-point vector stream function correlation has been devel- 
oped using tensor invariant arguments and evaluated by the comparison of model 
predictions with DNS data for incompressible homogeneous turbulent shear flow. 
This two-point vector stream function model correlation can then be used to cal- 
culate the two-point velocity correlation function and other quantities useful in 
turbulence modeling. The model assumes that the two-point vector stream func- 
tion correlation can be written in terms of the separation vector and a new tensor 
function that depends only on the magnitude of the separation vector. The model 
has a single free model coefficient, which has been chosen by comparison with the 
DNS data. The relative error of the model predictions of the two-point vector 
stream function correlation is only a few percent for a broad range of the model 
coefficient. Predictions of the derivatives of this correlation, which sure of interest 
in turbulence modeling, may not be this accurate. 


1. Introduction 

In one-point modeling for second-moment closure, four different terms are un- 
known: the pressure-strain term the dissipation term £,*j, the pressure diffusion 
term p Dij , and the turbulent diffusion term l Dij. A large amount of information 
that is needed for closure is contained in the two-point correlation tensor Rij. Given 
a model for i?ij, one can express the Reynolds stress tensor, the dissipation tensor, 
and the rapid part of the pressure-strain tensor as functions of Rij: 


u\u) = Rij(x,r = 0) , 


£a = lim v 
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where (i <-► j) in Eq. (3) indicates the addition of the previous term with inter- 
changed indices i and j. Rij also provides turbulence length-scale information, 
including the integral length scale and the Taylor microscale. 
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For engineering applications it is impractical to solve the two-point correlation 
equation to get the required one-point information. The approach should be to 
extract information from the two-point correlation tensor and the two-point cor- 
relation equation in order to improve one-point models. A first step was provided 
by the two-point correlation equation for isotropic turbulence developed by von 
Karman &: Howarth (1939), which gave early insight into the decay of isotropic 
turbulence and the development of turbulence length scales. This isotropic tensor 
form was used by Crow (1969) to derive an exact expression for the rapid pressure- 
strain- rate correlation in isotropic turbulence. A crucial further step was taken by 
Naot et al. (1973), who extended the two-point correlation model for i?,j, which 
led to the linear rapid pressure- strain model utilized in the well known Launder, 
Reece and Rodi (1975) second-moment-closure model. 

One-point turbulence models require some sort of length-scale equation. The first 
complete Reynolds stress model was developed by Rotta (195 la, b) in combination 
with a scalar integral length-scale equation. He developed the length-scale transport 
equation by introducing an integral operator to the trace of the two-point correlation 
equation. Wolfshtein (1970) has developed a similar equation. Neither Rotta or 
Wolfshtein introduced a model for the two-point correlation Rij, but they did model 
all the sink, source, and diffusion terms on the right hand side of the length-scale 
equation empirically. 

A completely different procedure was introduced by Donaldson and co-workers 
(Sandri (1977, 1978), Sandri & Cerasoli (1981), Donaldson & Sandri (1981)), who 
developed a new tensor length-scale equation. They approximated the correlation 
function R,j as a delta peak at zero separation and applied an integral operator 
to the transport equation for Rij . This crude assumption leads to the elimination 
of some important terms in the tensor length-scale equation. In addition, they 
introduced models for some terms where exact terms can be derived. 

Recently Oberlack ( 1994a, b) developed a new tensor length-scale equation. In 
his approach he has introduced a new model for the two-point correlation equation 
on the basis of tensor invariant theory. With this new approach the linear part of 
the rapid pressure-strain model could also be rederived. 

The objective of this paper is twofold. First, the model introduced by Oberlack 
(1994a) in terms of a tensor potential will be recast in terms of the two- point 
correlation of the vector stream function and generalized to facilitate comparison 
with the DNS results. Second, the basic model assumption made by Oberlack 
will be evaluated by comparing the model for the two-point vector stream function 
correlation with DNS data of a homogeneous shear flow calculated from Rogers et 
al. (1986). 

Section 2 gives an introduction to the vector stream function and discusses the 
non-unique behavior resulting from its definition, and then gives the definition of the 
basic two-point variables in correlation space and the relationship between them. 
Section 3 focuses on the derivation of the two-point vector stream function correla- 
tion model. A procedure to compare the model with the DNS data is outlined, and 
the reduction to a one parameter family of model equations is explained. Section 4 
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presents the results. There we introduce a new scalar measure that quantifies the 
difference between the DNS data and the model predictions. This quantity is used 
in determining the single model parameter 7. 

2. Definition of the vector stream function 

The new model for the two-point velocity correlation function introduced by 
Oberlack (1994a, b) had to be redefined for comparison with the DNS data of Rogers 
et al (1986). For this purpose it is helpful to start with the vector stream function 
formulation. 

The turbulent fluctuation velocity u\ can be expressed as the curl of the vector 
stream function fluctuation (Aris 1962), 


, dil>' k 

Ui = eiik di~ ’ 


(4) 


where f tJ k is the alternating tensor. 

Obviously any given determines u\ uniquely, but not vice versa because any 
curl- free vector field t/>* can be added to V>' without changing the velocity vector u\. 
Thus, Eq. (4) does not define a unique ip [ , and an additional condition is required. 
Talcing the curl of Eq. (4), we get the Poisson-like equation 


, _ d 2 ip' k 

U> dxidxk dxkdxk ’ 

where is the turbulent vorticity fluctuation given by 


(5) 


Wj = tijk 


du[ 
dxj 


( 6 ) 


Eq. (5) is fundamentally different from the standard Poisson equation, a fact that 
is easily seen in wavespace where the derivatives in (5) are replaced by multiplica- 
tion by wavenumber components. In wavespace, the Poisson equation is a regular 
linear algebraic equation with a unique solution for i/>J, whereas (5) is a singular 
linear algebraic system that requires the specification of an additional constraint 
to determine t/>J. To define a well- posed problem, it is necessary to introduce an 
additional restriction on as noted above. 

A numerically very useful and often used condition on V’l is that it be solenoidal 


M 

dxi 


= 0 . 


By using (7), Eq. (5) reduces to the usual Poisson equation 


(7) 



ay; 

dxkdxk 


( 8 ) 


Eq. (8) has been used to determine the vector stream function from the DNS 
data. However, this implicit use of (7) makes it necessary to reformulate the model 
of the tensor potential given by Oberlack (1994a). 
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All the two-point correlations functions below follow the definitions of Oberlack 
(1994a, b). We introduce the spatial vectors x and f* 1 * and define the correlation- 
space vector as 


f = i* 1 ) — x . (9) 

In general, the two-point correlation functions depend on the physical and the 
correlation space coordinates x and r and on the time t. 

The tensor potential V mn introduced by Oberlack (1994a) can be expressed in 
terms of the two-point correlation of the vector stream function fluctuation 


V mn (x, r, t ) = t/>' m {x, ^{.(xW, t ) . (10) 

V mn contains important length-scale information of the turbulence and can be di- 
rectly linked to the two-point velocity correlation function 


R i} (x,r,t) = u' i (x,t)u' j (x^ 1 ') ,t) . (11) 

R,j is the basic quantity in two-point modeling. The link between V mn and R tJ can 
be found using Eq. (4) and the definition of the correlation space r given in Eq. (9), 
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€ikm Cjlrt 
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(12) 


The tensor potential V mn was originally introduced to form an R{ } that would 
automaticEilly satisfy the two equations 


dRjj 

dii 


= 0 and ^>=0 


dr{ 


dr } 


(13) 


emerging from the continuity condition for the turbulent velocity fluctuation 
The constraint for the vector stream function %l>\ in Eq. (7) imposes two additional 
restrictions on the two-point vector stream function correlation 


dVg 

dxi 


dV H n . dV tj 

= 0 and — 0 . 

vTj dv j 


(14) 


Oberlack (1994a, b) did not express V mn in terms of and hence the constraints 
(14) were not implemented. It can be shown that he used an implicit constraint on 
the two-point vector stream function, but this is not computationally useful. 

An important consequence of introducing (7) is that it becomes necessary to 
reformulate the model for the two-point vector stream function correlation because 
the model of Oberlack is not able to satisfy (14). The derivation of the generalized 
model is described in the next section. 
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3. Derivation of the two-point vector stream function correlation model 

The ideas used here to model the two-point vector stream function correlation 
are similar to those introduced by Oberlack (1994a), but additional terms have 
been introduced to generalize the formulation. Here only homogeneous flows in in- 
finite domains are considered, hence statistical quantities are independent of spatial 
location and all the spatial derivatives in Eqs. (12), (13), and (14) vanish. 

We seek to model the two-point vector stream function correlation V mn in terms 
of a new tensor function G mn and the separation vector r. The components of G m n 
are all assumed to be functions only of the magnitude of the separation distance 
r = |r| and time t, 


G m „ = G mn (r,t) . (15) 

Because of the assumption (15), the tensor G mn has no angular variation in corre- 
lation space; the value of each component is constant on spherical shells in r-space. 
For simplicity, G mn has also been assumed to be a symmetric tensor. Thus, for 
general homogeneous flows, G mn is characterized by six independent functions of r 
and t . This is in contrast to isotropic turbulence, which is characterized by a single 
scalar function f(r y t) (see Appendix). 

The most general V mn that is linear and non-differential in the tensor G mn will 
involve the following terms: 


^ r ^ c nbn r mTn r n r k r m r n r k r, 

Grnm vmnG kk > Vmn ^ Gkh ^2 G kk y ^2 ^ 2 ^ 4 kl • \ ) 

However, (14) can not be satisfied by an arbitrary linear combination of these forms. 
Therefore, we add an additional set where G m n I s replaced by 
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^ mn * ’ q 
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The resulting generalized form of the model is then 
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The coefficients a, and 0, are assumed to be constants for any given field. The 
tensor Vm° n contains the value of V mn at r = 0 and has to be introduced to avoid a 
singularity at the origin. It does not contribute to R,j or other quantities of interest 

derived from the model because jffrVmn = 0. Throughout the rest of this paper, 

we have absorbed Vmn into V mn so all elements of V mn drop to zero at r = 0 (t.e., 
we consider the quantities V mn - V^). 

From its definition, V mn is Hermitian in homogeneous turbulence; that is V m „(r) 
= u — r). The model Eq. (18) is more limited and predicts that V mn is symmetric 
since G mn has been assumed to be symmetric. 

Because the tensors V mn and G mn are symmetric, only one condition in (14) has to 
be satisfied. Substituting the model (18) into (14) results in an equation containing 
nine independent vectors. The nine coefficients of these vectors are linear functions 
of a, and j3 % and are set to zero to satisfy (14). This ensures that the model two- 
point vector stream function correlation V mn is consistent with the assumption of 
a solenoidal vector stream function. 

We thus have 9 equations for the 12 coefficients a* and Without any loss of 
generality, we have chosen as, /? 4 , and /3$ to be the independent coefficients and have 
solved the equations to express the remaining nine coefficients as linear functions 
of these three parameters. The solution yields 


11 3 1 

a i = ~2 a5 ~ > Q2 = 2 a * “ 2 ^ 4 ” ^ 5 ’ a3 = - 2 as ’ “ 4 = 

Q 6 — — 2 Q 5, 01 = ~05, 02 = —04, 03 = — ~ 05, and 06 = ^<*5 . (19) 

At this point, the DNS data of Rogers et al. (1986) can be used to optimize the 
choice of the model constants a 5 , /? 4 , and /9 5 . From the DNS data, V mn can be 
calculated directly, but we also need a procedure for evaluating the six independent 
functions G mn in order to assess the quality of the model. 

Eq. (18) consists of a coupled set of ordinary differential equations for G mn . (Note 
that in previous formulations of the model that did not have the r terms, 

G mn could be obtained by solving a system of six linear algebraic equations.) The 
model does not specify the functional dependence of G mn on r, and the DNS data 
must be used to evaluate these functions. 

We have done this by equating the spherical- shell- averaged V mn , computed from 
the DNS data, with the spherical-shell-averaged V mn predicted by the model using 
(18). Applying the definition of spherical-shell averaging given by 

2w n 

v mn = V mn dV = ^jj v m „ sm(6)d8d<f> 
n oo 

to the model expression in Eq. (18) with the constants defined in (19) yields 


( 20 ) 
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V’mn — (30 tt5 3 ^ 5 )( 3G! mn+r ^ ) 

+ (~l 5 a5 " l 13 * ~ 3 A ) * m " ( 3G “ + r_ aT ) - (21) 

This coupled system of ordinary differential equations can be solved by splitting 
G mn into its trace Gkk and its anisotropic part g mn — G mn /Gkk — ^mn/ 3 . This de- 
composition decouples the differential equations for Gkk and those for each element 
of gmn' Solution of these uncoupled equations then results in the determination of 
G mn (r) in terms of the spherical-shell-averaged V mn 

t r 

G mn = J V mn r”dr' + 4 ^ J V k kr' 2 dr' , (22) 

0 0 

where 

Al “ a 5 - 1 O 0 5 ^ Ai a 5 + 12& + 80 s a 5 - 1O0 5 ' ( 1 

All the integration constants have been set to zero to eliminate a singularity at 
r = 0. 

For a given set of 05, /?4, and /J5, the functions G mn can be computed from the 
spherical- shell- averaged V mn obtained from the DNS data. At this point it is clear 
that the functions G mn computed in this manner do not depend on all of 05, /?4, 
and ft, but only on the reduced parameter set (Ai,^). 

To distinguish clearly between V mn predicted by the model (using the functions 
Gmn computed from the DNS data as described above) and V mn calculated from the 
DNS data, superscripts are used. The vector stream function two-point correlation 
calculated from the DNS data will be labeled as and that predicted by the 

model will be termed Substituting G mn as determined from (22) into the 

model ( 18 ) and ( 19 ) results in the following model expression: 


T/ r model 
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where 

r 

Rmn = pj- J V mn r' 2 dr' and 7 = . (25) 

0 

Note that the model (24) does not now depend on q 4 and that q 5 and /? 5 only 
appear in the ratio 7 . Thus the model has only one free parameter that remains to 
be chosen by comparison with DNS data! 

In the next section we compare the values of V mn computed directly from the 
DNS with those predicted by the model given by (24). This comparison is made 
for many values of 7 and permits the determination of an optimum value of 7 that 
results in the best agreement between the model predictions and the DNS data. 


4. Results 

As noted above, direct numerical simulations of incompressible homogeneous tur- 
bulent shear flow generated by Rogers et al (1986) have been used both to optimize 
the choice of the model constant 7 and to test the validity of the model assumptions 
by comparing model predictions_with DNS results. In these numerical simulations, 
the streamwise mean velocity U varies linearly in the cross-stream (y) direction, 
with the constant and uniform mean shear rate being given by S = dU jdy. In 
this work the C128U simulation has been used. This case evolves until a non- 
dimensional time of St = 16, beyond which the computational domain becomes too 
small to capture an adequate sample of large-scale energetic turbulent eddies. 

Although fields at various times St have been examined, results presented here are 
for St = 10. Results at other times are qualitatively similar, but after this time the 
stream function (primarily from eddies of larger scale than the energy- containing 
eddies) begins to be limited by the computational domain size. The flow cannot be 
regarded as a developed shear flow before about St = 8 . 

As noted earlier, the model has been constrained to predict symmetric vector 
stream function two-point correlations V mn . Therefore the DNS results have been 
symmetrized prior to being used for model development and evaluation (this results 
in only minor adjustments to the off-diagonal tensor components of V mn ). 

The degree to which V^ s and V™° dcl are different is a measure of the quality of 
the model. A simple measure to quantify this difference is given by 


„ Jy<* 

p = 


- V™ 5 ) 2 <Pr 


f v (V™ s f<P r 


(26) 


where the integration domain V is restricted to a sphere of radius r mar , the largest 
radius that will fit into the computational domain (r maz = min(L x /2, L y j 2 , I*/ 2 ), 
where L, is the computational domain size (periodicity length) in the i th direction). 
Integration over larger separation distances is not possible because the data would 
be contaminated by periodic images in the computation. Note that each point in 
correlation space is weighted equally in the above definition. Alternative measures 
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FIGURE 1. Dependence of the relative error p on the parameter 7 for a) positive 
values of 7 and b) negative values of 7. Note that there is a singularity at 7 = 0.1. 

that weight the data for smaller separation distances more heavily do not change 
the optimum value of the model constant 7 significantly. 

Fig. 1 illustrates the dependence of p on the model constant 7 for both positive 
and negative values of 7. According to this measure, the model is equally good for 
both large positive and large negative values of 7, with a value of about 0.018. There 
is a singularity at 7 = 0.1. Of the values of 7 examined (indicated by solid circles 
in the figure), the lowest value of p = 0.0172 is achieved at 7 = -20. However, 
there is a broad minimum with many other values of 7 yielding values of p that are 
almost this low. Other considerations discussed below suggest an optimum value 
of 7 =3 —5, with p = 0.021. Figs. 2 and 4 have been generated using this value of 
7 = -5- 

The qualitative form of p(7) for the fields at St = 8 and St = 12 is similar to 
that shown in Fig. 1, although the value of p at the broad minimum is different. At 
St = 8, there is a minimum of p — 0.0055 at 7 = —20; by St — 12, p increases to 
0.060 at 7 = —20. By St = 12 the computational domain is no longer sufficiently 
large to capture V33 and the increased value of p is therefore not surprising. 

The functions G m „, derived as outlined in §3, are shown in Fig. 2 for the C128U 
simulation at St = 10 using 7 = —5 (Ai = 20/17, A 2 = —800/969). Because Vm„ 
has been subtracted from the vector stream function two-point correlation to obtain 
V mn , V m n starts at zero for r = 0 and approaches a constant value of — Vm°n for 
large separation distances where the correlation falls to zero. Using Eq. (22), it can 
be shown from this that G mn is also zero for r = 0 and that the functions G m n will 
also approach constant values as the separation distance r becomes large. From 
Fig. 2 it is clear that the tensor components G 12 and G 22 level off by r — r mai , 
whereas G u and G33 do not. The off-diagonal components G13 and <? 2 3 are nearly 
zero compared to the other components (similar behavior is observed in the one- 
point correlations u\ u' :i and u' 2 u' 3 , which are also zero in homogeneous turbulent 
shear flow). 

The measure p indicates what good choices for the model constant 7 are, but it 
is also instructive to examine directly, and to compare it with Vf 2 NS . Fig. 3 

contains contour plots of the six independent V££ s components. All contour plots 
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FIGURE 2. Functional dependence of G mn on r for the C128U simulation at 
St = 10 and 7 = —5 (A x = 20/17, A 2 = -800/969). 


shown are in r x -r y -planes at r z = 0. As a result of the mean shear, the contours 
are elongated in preferred directions. The “ridges” of the V^* 5 , V^ NS ^ and V™ 3 
contours are inclined to the streamwise direction at an angle of somewhat less than 
45°, as is often the case with many quantities in this flow. The V^ NS contours are 
qualitatively different, with the “ridge” running roughly orthogonal to that of the 
other components. Similar behavior is observed at other times St, The off-diagonal 
terms Vi 3 * s and V^ NS are an order of magnitude smaller than the other terms and 
exhibit less organized structure; presumably they would be zero with an improved 
statistical sample of eddies. Examination of fields at earlier St lends credence to this 
suggestion. At earlier times the eddies are smaller, resulting in an increased sample 
and improved statistics. For these fields the V^ NS and V^ NS components are about 
two orders of magnitude smaller than the others and about the same magnitude 
as the level of the fluctuations for large r in the other components (unlike for the 
flow at St = 10, at earlier times the V£„ s contours are confined to a much smaller 
region in the center of the computational domain). 

Although the measure p is insensitive to the value of 7 in the range 7 < —10 
and 7 > 10, the orientation of the contour lines for V™° del varies significantly. The 
contours of for the C128U simulation at St = 10 for 7 = — 5 are illustrated 

in Fig. 4. The model predictions are limited to the interior of a sphere of radius 
Tmaz, because the functions G mn cannot be determined from the DNS data for 
separations r greater than this as noted above. For this choice of 7 , the contour 
ridge orientations of the various components are approximately the same as those 
of the V£” s shown in figure 3. The peak values and their decay in all directions are 
reasonably well represented within the radius considered. Beyond this radius, the 
DNS data show an increasing complexity of structures with additional secondary 
maxima that may be strongly influenced by the limited statistical sample of the 
computation and are not predicted by the model. 

The orientation of the contour ridges for other values of 7 is different. When 7 is 
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Figure 3. Iso-contours of a) b) c) V& NS , d) V 2 D 2 » S , e) V 2 ° NS , 

and f) V£ NS for the C128U simulation at St = 10 in the r^-r^-plane at r z = 0 . 
Value at the center of the domain ( r x = 0 and r y = 0) is zero, solid contours 
indicate negative contour levels, and dotted contours indicate positive contour levels. 
Contour increments are a) 0.004, b) 0.005, d) 0.010, f) 0.010, and contours levels 
in c) axe -0.006, - 0 . 002 , and 0.002, and in e) are -0.014, -0.010, -0.006, -0.002, and 
0 . 002 . 

decreased to larger negative values, the V{[° dcl and contours become aligned 

with the coordinate directions, V 1 7° d ' 1 being oriented in the streamwise direction 
and V 22 odel being oriented vertically. The V 12 ,,d “ and V£° d “ contours do not change 
significantly. 

The V™° d ‘' contours are the same for both large positive and large negative values 
of 7 , but for smaller positive 7 the orientation is different from that associated with 
negative 7 . Again, the V{ 2 ° del and V 3 ™" d ' 1 components are relatively unaffected by 
changes in 7 , but for 7 = 2 the contour ridges of the Vfi° d '‘ and V 2 Y d “ components 
are opposite those shown in Fig. 4, i.e. the V 7 1 7 ° Je ' contours are inclined at a slightly 
negative angle to the streamwise direction rather than positive, and the V 22 ode ' 
contours are inclined to the left of the vertical axis rather than to the right. For 
larger positive 7 the coordinate axis-aligned contours associated with large 7 are 
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Figure 4. Iso-contours of a) V{{ odcl , b) V{$ odct , c) d) Vfi odtl , e) 

and f) V^ odet for the C128U simulation at St = 10 in the r^-r^-plane at r z = 0 
for 7 = — 5. Value at the center of the domain ( r x = 0 and r y = 0) is zero, solid 
contours indicate negative contour levels, and dotted contours indicate positive 
contour levels. Contour levels are the same as in Fig. 3. 

rapidly approached. 

It is thus apparent that negative values of 7 should be chosen to best capture the 
contour orientation found in the DNS results. In this respect, this criterion leads 
to a similar conclusion as that based on the measure p (which reaches a minimum 
for negative 7). However, while the measure p suggests an optimum value of about 
7 = —20, matching the contour orientation would lead to a different choice of 7 
(smaller in magnitude). From Eq. (12) it is clear that to predict accurately the 
two- point velocity correlation function (and thus the Reynolds stresses) we need to 
be able to obtain accurate estimates of derivatives of V mn . Simple agreement in 
the magnitude of V mn , as measured by p, does not ensure good prediction of the 
Reynolds stresses. Since virtually all quantities of interest to be computed from 
the model require differentiation of V™* del to obtain, we felt that sacrificing a few 
tenths of a percent in p in order to get the contour orientation correctly would be 
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justified and have taken 7 = -5 as the optimum model constant as noted above. 
Ideally one would like to optimize the choice of 7 by comparing model predictions 
for the quantity of interest ( e.g . Reynolds stress) directly with the DNS values of 
the same quantity. It is possible that in other flows the optimum value of 7 would 
be different. 

5. Conclusions 

A new model for the two-point vector stream function correlation has been de- 
veloped. The proposed model has been compared with DNS data for homogeneous 
turbulent shear flow. The model gives a relative error p of a few percent for a broad 
range of the single model parameter 7. A more limited range of 7 also yields the 
correct orientation of the two-point vector stream function correlation contours. 
Any integral variables (like the integral length scale) calculated from the model 
Eire relatively insensitive to the parameter 7 and should model the DNS data quite 
well. The model was developed for the purpose of modeling the two-point velocity 
correlation. 

The results presented here suggest that the fundamental model assumption, 
namely that the two-point vector stream function correlation can be modeled in 
terms of the separation vector and a symmetric tensor Gmn whose components are 
only a function of the magnitude of the separation vector, is justified. However, if 
more directional information is needed and differentiation is applied to the model 
predictions, the resulting quantities may have some shortcomings. An example is 
the rapid part of the pressure-strain correlation, which can be calculated from the 
model but leads to the usual linear model with its known deficiencies. 
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Appendix 

For homogeneous isotropic turbulence the tensor G mn can be written in terms 
of a single function G(r) as G mn = G{r)6 mn . Substituting this into Eq. (18) and 
using (19) yields 


1/(0 + — (Al) 

V mn - F mn - (r) + 2dr J 6mn r2 2 dr 1 ^ 


where F = -(a 5 /2 + 6/? 4 +4/3 5 )G. Note that this is of the same form as the standard 
expression for the two-point velocity correlation tensor in isotropic turbulence 


R,, = I* [(/(r) + r -f(rj) 6„ - 5/V)^] ■ (A2) 

Employing Eq. (12), we can relate the function F(r) in equation (Al) to the function 
/(r) in Eq. (A2) 
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F(r) = f \h[ r * fdr '- [ r ' fdr ' ■ 


(A3) 
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Rotating turbulent flows and modeling 

The motivation for the group effort in this section stems from the fact that it 
can be shown that all one point closure models currently used in RANS (Reynolds- 
averaged Navier-Stokes) codes fail to correctly model the effects of strong rotation 
on turbulence. Yet, turbulent flows that are subjected to system rotation are found 
in many situations, examples are turbomachines, the wing tip vortex, geophysical 
flows, and many others. A prime objective of the CTR is to foster a climate where 
new ideas on understanding and modeling of turbulent flows will generate innova- 
tions that become part of the foundation of the turbulence knowledge base. The 
five papers in this section form a good example where innovations in modeling and 
new insights into the effects of rotation on turbulence are achieved through the use 
of a combination of both theory and simulations (direct and large-eddy). 

Blaisdell and Shariff performed direct numerical simulations of homogeneous tur- 
bulence in elliptic streamline flow. This is a new flow in a series of building block 
flows where particular issues about turbulence modeling and physics are addressed. 
In this flow, the effects of both rotation and strain combine to provide an instability 
mechanism that has been proposed as a universal mechanism for energy transfer 
from large scales to small scales. Four cases were simulated to investigate the effects 
of Rossby number and ellipticity. A fifth case at high resolution was carried out 
to test resolution and sampling issues. Statistics of interest to turbulence modeling 
are presented. The simulations lead to the discovery of a homogeneous flow where 
the non-linear cascade is periodically suppressed and re-established. 

Mahalov derived and analyzed the evolution equations for long-time averaged ro- 
tating shallow- water equations. He then used the same approach to derive equations 
for homogeneous flows subjected to background rotation. In this case, he decom- 
posed the flow field into two-dimensional modes that are unaffected by rotation and 
three-dimensional disturbances. He showed a connection between his physical space 
formalism and the helical wave space formalism. The derived equations may prove 
useful in explaining the effects of rotation on the turbulent transfer between the 
scales. 

Squires, Chasnov and Mansour used large-eddy simulations to investigate the 
asymptotic similarity of rotating homogeneous turbulence. They build on their pre- 
vious investigations where, in the limit of high Reynolds number and small Rossby 
number, power laws for the turbulent stresses and length scales were derived based 
on dimensional analysis and simulation results. In the present study, a search for 
similarity laws for the spectra was carried out. They defined four independent en- 
ergy spectra and found scalings that would collapse three out of the four spectra. 
The spectrum which does not collapse corresponds to two component motions in 
the plane normal to the rotation axis. A strong reverse cascade of the energy from 
small- to-large scales is found. These results hint at the tendency of the flow to 
become two-dimensional. 
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Cambon, Mansour and Squires used the same large-eddy simulations as above 
to investigate the development of anisotropies and the tendency of the flow to 
become two-dimensional. This tendency is often used as a confirmation of the 
Taylor- Proudman theorem. However, it can be shown that a key assumption in the 
derivation of the theorem does not hold for homogeneous flows. Cambon has long 
argued that the observed two-dimensionalization in rotating flows is a non-linear 
phenomenon. In the present work, results from low- Reynolds number direct numer- 
ical simulations and from high- Reynolds number large-eddy simulation are unified 
by arguing for the existence of two transitions in rotating turbulence that can be 
identified by two key Rossby numbers. During the first transition, anisotropies 
develop when a macro Rossby number drops below unity. During the second tran- 
sition, the nonlinear transfer is shut off when a micro Rossby number drops below 
unity. The large-eddy simulation fields were also analyzed for alignment between 
vorticity fluctuations and the rotation axis. The results confirm the dominance of 
corrotative vorticity in the long-time high-Reynolds number limit. 

Hadid, Mansour and Zeman tested a new one-point closure model that incor- 
porates the effects of rotation on the power-law decay exponent of the turbulent 
kinetic energy. This modification to the e-equation was proposed by Zeman based 
on the large-eddy simulation results. During the summer program the model was 
successfully tested against experimental data of isotropic turbulence that was sub- 
jected to uniform rotation. A new definition of the mean rotation was proposed 
based on critical point theory to generalize the effects of rotation on turbulence to 
arbitrary mean deformations. The model was then incorporated in a code used at 
Rocketdyne and tested for a backward-facing step and a dump combustor. Prelim- 
inary results show promise. Further collaborations on modeling rotating flows will 
continue between CTR, NASA/Ames, and Rocketdyne. 

N, N. Mansour 



